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Abstract 

We investigate a periodic version of the Benjamin-Ono (BO) equation as- 
sociated with a discrete Laplacian. We find some special solutions to this 
equation, and calculate the values of the first two integrals of motion Ii and 
I 2 corresponding to these solutions. It is found that there exists a strong 
resemblance between them and the spectra for the Macdonald g-difference 
operators. To better understand the connection between these classical and 
quantum integrable systems, we consider the special degenerate case corre- 
sponding to q = in more detail. Namely, we give general solutions to this 
degenerate periodic BO, obtain explicit formulas representing all the integrals 
of motions I n (n = 1, 2, • • •), and successfully identify it with the eigenvalues 
of Macdonald operators in the limit g — *> 0, i.e. the limit where Macdonald 
polynomials tend to the Hall-Littlewood polynomials. 

Key words: Benjamin-Ono equation, nonlocal integrable system, 
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1. Introduction 

In [3], the authors studied a doubly periodic version of the intermediate 
long wave (ILW) equation associated with a discrete Laplacian. Then this 
ILW-type equation was identified with the first member of an integrable 
hierarchy which is given as a certain reduction of the discrete KP theory in 
the framework of the Sato theory. 

In this paper, we study a spacial case having a single periodicity, instead of 
dealing with the most general doubly periodic case which requires much more 
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technique based on the algebro-geometric argument. Namely we consider the 
periodic Benjamin-Ono equation [2], 0] with discrete Laplacian 

d f 1/2 
—r](x,t)=r)(x,t)J ^ [cot{7r(y-x-7)} 

idy 

-2 cot {n(y -x)} + cot {n(y - x + 7)}] rj(y, t)— , (1) 

where the symbol f denotes the Cauchy principal value integral and 7 is 
a complex constant with nonzero imaginary part. Our goal in this paper is 
to study some special solutions to ([!]) and the properties of the associated 
integrals of motion in some detail. We also investigate the special case 7 — > 
zoo in great detail since it allows us to have very much explicit results about 
solutions and integrals. 

We claim that (pQ) can be regarded as a classical integrable system as- 
sociated with the Macdonald theory of symmetric functions with two pa- 
rameters q and t |4|. Recall that Macdonald introduced a set of commuting 
g-difference operators D\, D2, ■ ■ ■ depending on q and t, which are acting on 
the space of symmetric polynomials, say in x^'s. Then he proved the funda- 
mental existence theorem for the simultaneous eigenfunctions P\(x;q,t) of 
the DiS, together with a certain normalization condition. Here the index 
A = (Ai, A2, • • •) denotes the partition. 

For our purpose to compare (JTJ) with the Macdonald theory, we better 
consider the infinitely many variable case and use the Heisenberg represen- 
tation of the Macdonald operators developed in [AMOS], [S] and [FHHSY]. 
Introduce the Heisenberg algebra given by the generators {a n }nez^ an< ^ the 
commutation relations 

1 _ q\m\ 

[a m , a n ) = m- -— -5 m+n)0 . (2) 

1 — v ' 

Let £ be a constant. Set 

T){z) = ^2 VnZ~ n = e exp | ^ - — - — a_ n z n | exp | - ^ - — —a n z~ n | . (3) 

neZ \n>0 " /\ n>0 n J 

It was shown in [S] and [FHHSY] that we have an infinite commutative family 
of operators acting on the Fock space containing t]q as the first nontrivial 
member. 
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One of the canonical ways to write down the generators of this commu- 
tative algebra is as follows. Set 

/»=/• J] ^^■■■V(zi)---V(zn):, (4) 

J Jc 2mzi 2mz n *-\ Zj - qz k 

l<j<fc<n J 

where the integration contour C is the torus \zj\ = l,j = l,---,n, and 
the symbol : ■ • ■ : denotes the ordinary normal ordering of the Heisenberg 
generators, more explicitly we have 

l<j<fc<« 1 3 kA3 q k> 

Then we have [I m , I n ] = for all m, n > 1. Moreover we can write down ex- 
plicitly the spectrum of J„ as g-™^" 1 )/ 2 Y[ n k=l (l-q k )-e n (et- x \ Eqt~ X2 , eq 2 t- X:i , ■ ■ 
where e n (xi,X2, ■ ■ ■) denotes the n-th elementary symmetric function, A = 
(Ai, A 2 , ■ ■ ■) is a partition. 

For simplicity, we set a n = —(1 — t n )a n /n. With this notation, we have 

7]{z) = e : exp (^n^o a ™ 2;_n ) ■• Now we proceed to considering a classical 

limit. Set t = e h . While fixing q, we consider the limit h — > 0, namely t — *■ 1. 
In this limit we have [a m , a n ] = and [r) m ,r) n ] = 0. Hence we regard the 
algebra generated by a n 's or r/^s being our algebra of classical observables. 

Induce the poisson bracket by {u, v} := limf^o[u,v]/H as usual. Then 
we have the set of canonical commutation relations {a m ,a n } = sgn(m)(l — 
g H )5 m +n,o, which gives us {r] m ,Vn} = £ 1? to sgn(Z) (1 - q W )v 

Let ?7o be our Hamiltonian. For example, the time evolution of r](z) is 
given by 

d 

— rj(z) = {rjo, rj(z)} = r)(z) sgn(n)(l - q lnl )r]„ n z n , (6) 

which is just identical to (pQ) under the identification of the variables z = 
exp(27rix) and q = exp(2mj). 
Set 

t+(z) = exp I - Y^ a ~ nZn ) ' = 6XP ( ~ S iL ""^ ) ' ( ? ) 

V n>0 ^ / \ n>0 ^ / 
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We have 




(8) 



D t T-{z)- t + {z) = er_(g l z)r + {qz) - rj T_{z)T + {z) . 



(9) 



Here comes a question: how one can find some special solutions to ([T]) or 
([9]), calculate integrals of motion, and find some connection with the theory 
of Macdonald symmetric polynomials? In this paper, we show some explicit 
connection between the (classical) integrals and the (quantum) eigenvalues, 
indicating that some deep structure is hidden in our problem. 

This paper is organized as follows. In §2, we give some special solutions 
to ([T]). The first few J n 's of them are calculated explicitly. It is shown that 
they have precisely the same formulas as the ones for the eigenvalues of the 
Macdonald operators. In §3, the special case q = is treated in a complete 
manner. The general formula for J n 's is given. 

A remark is in order here. In the work by Avanov, Bettelheim and 
Wiegmann Q, a classical integrable system associated with the Calogero- 
Sutherland model is studied. It is an interesting future problem to connect 
their bilinear equation and solutions with ours. 

2. Some special solutions to the periodic BO equation (TTJ) 

After a little algebra, one can construct some Laurent polynomial solu- 
tions to the bilinear equation ([9]). Recall that e is a parameter in (Q, and 
we will consider it as an arbitrary parameter. We have a class of solutions 
parametrized by one more parameter a. One may easily find r + and r_ given 



satisfy (jUJ) under the condition that eq + (1 — q)a — r)o — I\. 

We can compute higher integrals of motion (j4]) corresponding to this 
special solution. For example, we have 



by 




(10) 



h = eq + (1 - q)a = (1 - q)e 1 (a, eq, eq 2 , ■■■), 




e 2 (a,eq,eq 2 , ■ ■ ■), 



(11) 



4 



and so on. 

Next, we show another class of special solutions to 
parameters a\ and a 2 . Set 



which has two 



1 + ze {1 - q)ait + ze (1 - q)a2t 



+ 



(ai - a 2 f 



(ai - ga 2 )(ai - q 1 a 2 ) 
r_ = 1 + c 1 z- 1 e-( 1 - q)ait 
4 



z 2 e (l-g)(ai+a 2 )t 



C2Z -l e -(l-q)a 2 t 



;i2) 



(ai - a 2 



(ai - ga2)(a-i - q 1 a 2 ) 
e - qaj (ai - qa 2 )(a 1 - g _1 a 2 ) 



Cl c 2 z- 2 e~ {1 - q){ai+a2)t , 



e — q z cij 



(a : - a 2 ) 2 



(13) 



Then the ([9]) holds under the condition that eq 2 + (1 — q)ai + (1 — q)a 2 = 
1] = Ii. One can calculate the higher integrals of motion The first two 
read 



h = eq 2 + (1 - q)a x + (1 - g)a 2 = (1 - q)e 1 (a 1 , a 2 , eq 2 , eq 3 



[l-q)(l-q 2 
Q 



2„2 



-aia 2 + eq(l - q)a\ + eq(l - q)a 2 + e q 



(14) 



[l-q)(l-q 2 ) , 2 _ fl 3 

— e<2\a\, a 2 , eq ,eq , 



q 



and so on. 

It is clearly seen from these examples that there is a strong resemblance 
between such formulas for the integrals of motion and the eigenvalues of the 
Macdonald operators q ~ n ( n ~ 1 )/ 2 ]\ n k=l (l - q k ) ■ e n (et~ x \ eqt~ x \ eq 2 t~ x \ ■ ■ •). 
Namely, if we take the limit t — > 1, Aj — > 00 in such a way that we have the 
finite limits limeg i_1 t _Al = aj, then we recover the above formulas for the 
integrals of motion. At present the reason of this beautiful correspondence 
has not been investigated. 



3. the spacial case q = 

At present, unfortunately, our study on the special solutions to Qj and 
whose integrals of motion still remains primitive and heuristic. In the spe- 
cial limit q = 0, however, we can easily complete our program within the 
technique of linear algebra as we will show below. 
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Before we embark on the study of the bilinear equation, let us try to solve 
naively in the limit q = 0, namely 

d 

t^v( z ) = V( z ) s & n ( n )V- n z n ■ (15) 

Set 7]±(z) := Xln>o VTn z±n - Then we have the split equations 
d d 

-j.r]+{z) = r)+(z)(r) + (z) + r) ), ^:V-( Z ) = -V-( Z )(V-( Z ) + Vo)- (16) 
One immediately finds the general solution, 

f , _ -rftc + (z)e*>' ( , _ -yoc4z)e-^ 

V+[Z) d + (z) + c+(z)e^ V - {Z) d4z)+c4z)e-^ [U ) 

where d±(z), c±(z) are arbitrary holomorphic functions in some neighbor- 
hoods of z ±l = 0. 

One may ask if d±(z),c±(z) can be restricted to Laurent polynomials in 
z ±x , and the denominators d±(z) + c±(z) exp(±r] t) play the role of the tau 
functions t±(z). The answer indeed is yes. 

In the case q = 0, we have 

V(z)=£ ( ! , v (18) 
A t-(z) -T+{z)=e- rj T„(z)T+{z). (19) 

Setting t±(z) = d±(z) + c±(z) exp(±r] t) , one finds that (|T9|) reduces to the 
functional equation 

d+(z)d-(z) - c+(z)c-(z) = e/rjo. (20) 

Proposition. Let n be a nonnegative integer, and let e±i , • • • , e± n be 2n 
parameters. Let e m (z,t) = e m z m exp{sgn(m)r] t} for m <G Z^ , and define 
t±(z, t; e±i, - • • , e±n) = r±(«) by 

r+0) = 1+ Y e i+ e i2 e_ h + ^ e i3 e_ i2 e il + 

l<i<n l<ii<i'2<n I<ii<i2<i3<n 

H h e n e_ n+ ie n _2 • • -e_(_i)n, (21) 

r_(2) = 1 + e_j + e-^e^ + e- i3 e i2 e- h + 

l<i<n l<ii<i2<n I<ii<i2<i3<n 

H h e_ n e n _ie_ n+2 ■ • -e(_i)». (22) 
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Then these satisfy (TT9~I) under the condition that n^=i(l ~ e i e -i) — e /Vo- 
Proof. We have 



d+(z) = 1 + 2j e i 2 e -n^ 

l<il<i2<n 



C4- u 



E ^ + E 



f e ■ e ■ z iz ~ i2+il 4- • • • 



l<i<n 

d-(z) = i+ Yl 



I<il<i2<i3<i 

-22+21 



+ 



l<il<i2<n 



C-(Z 



E 



l<j<n 

and need to show that 



+ E 

I<ii<i2<i3<n 



d + (z)d_(z) - c + (z)c^(z) = - e,e_ 



A way to write d± and c± is 

d+(z) c+(z) 
c-(z) d-(z) 



i=i 



1 e n z T 
e- n z- n 1 



1 €\Z 



e^z 



The determinant of the above gives the relation (T27I) . 



(23) 
(24) 
(25) 
(26) 



(27) 



(28) 



□ 



Remark. For n > 0, and e±i, • ■ ■ ,e± n , we have constructed a solution to 
ffT9]) and have denoted them by t±(z, t; e±i, ■ ■ ■ , e± n ). It is desirable, however, 
for our later purpose to have more flexible notation in which we can treat t± 
with different numbers of parameters within a unified notation. It is clear 
that from infinite sequences (e±i, e±2, ■ ■ •) with only finitely many nonzero 
parts, say e-y = for i > n, we can construct the tau functions as above. 
We will denote them by t±(z, t; e±±, €±2, ■ ■ •) without the need of referring to 
the V. 



When q = 0, the integral of motion (j4j) becomes the Toeplitz determinant 



Vo V-i 
Vi Vo 

7] n -i r) n -2 



V-n+1 
V-n+2 

Vo 



(29) 
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Theorem. Let r] be satisfying the condition n^i(l — e i e -i) — £ /Vo, an d 
let t±(z) be defined by (12"T1) and (I2"2"j) . Hence (TT^i) is satisfied. Then define rj n 
by f|T8|) or (fTTj) . Then we have 



4+i = % fc+1 (l-eie-i) fc (l-e 2 e 



sfc-l 



(1 - efce- 



(30) 



for /c > 1. 

Proof. We know that J ra 's are t independent. Thus we may set t = in 
what follows in order to make our computation simple. Set 



V± (z,0) 



-rj c±(z) 



d ± (z) + c ± (z) n>Q 
for notational simplicity. Then, our claim in Theorem is 



-1 

Z/X -1 



-1 



-l) fe+1 ri(l-^-/) 



\fc-Z+l 



z=l 



This can be proved by Lemma stated below as 



/ -1 U-l 
v x -1 



\ 



J 



( 1 o 



0\ 



rf 1 0--- 



(1 - eie_i) • • ■ (1 - e fc e_ fc ) 

-1 



Vrf 



1 / 



fjfc-i 



-1 



-(1 - eie_i) •••(!- e k e- k )Ik 



Lemma. For iV > 1, we have 



/ -1 V-l 
v x -1 



-(N-l) 



\ 



( 1 \ 



(31) 



(32) 



□ 



.(AO 



/ -(1 - eie_i) •••(!- e^e.Ar) \ 



y i/jv i/jv-i • • • 



8 



Here, rj" 1 -* is the coefficient of z ±j in t±(z, 0; e±±, e±2, • • • , e± m , 0, 0, • • •), namely 



t±(z, 0; e±i, e± 2 , • • • , e± m , 0, 0, • • •) = ^ , ±J 

j'=i 



(33) 



Sketch of the proof of Lemma. We have the following three properties, 
which can be proved by induction on N. 

1. u N is a polynomial depending only on e±i, e± 2 , • • • , e±(jv-i) and ejy, and 
v_ N is a polynomial depending only on e± 1 , e± 2 , • • • , e±(jv-i) and e_N- 



2. 



±7V 



3. We have 



= (1 - eie_i)(l - e 2 e_ 2 ) •••(!- e N -ie- N+1 ). 



( r<3 o 

T (N) (JV) 
'-JV+1 r -JV 



\ '-1 '-2 



o \ 



/ ^1 \ 

^2 



( « +1 \ / 



(JV) 
-Af+2 








where, is defined by 



d±(z, 0; e±i, • • • , e ±m , 0, 0, • • •) =: ^ cig 

j'=i 



(m V. 



□ 
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